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First  we  shall  state  the  problem  and  give  a  simple  method  for  its 
solution.  Then  we  shall  describe  the  empirical  method  used  to  relate 
the  solution  to  the  problem  of  blast  loading  an  air  earth  interface. 

In  circular  cylindrical  coordinates  (r,  0,  z)  a  homogeneous, 
Isotropic,  solid  medium  occupies  the  semi-infinite  region  0  £  z  <  oo  , 
and  a  homogeneous,  isotropic,  fluid  elastic  medium  occupies  the  seml- 
lnflnlte  region  -  oo  <  z  £  0.  An  Impulsive  point  source  of  pressure 
is  located  at  r  =  0,  z  »  -h. 

The  motion  in  the  solid  is  a  solution  to  the  linear  elastic  wave 
equation.  Integral  expressions  for  the  motion  have  been  derived  by 
Swing,  Jardetzky  and  Press  (Elastic  Waves  in  Layered  Media,  New  York, 
McGraw-Hill ,  1937),  Strick  ("Propagation  of  Elastic  Wave  Motion  from 
an  Impulsive  Source  along  a  Fluid/Solid  Interface",  Philosophical 
Transactions  of  the  Royal  Society,  Series  A,  Vol.  251,  1959)  and  others. 
Therefore,  no  derivation  of  these  expressions  will  be  given  here. 

Let  Q  be  the  r  component  of  particle  velocity,  W  be  the  z  component 
of  particle  velocity,  and  TIZ,  Tzr,  Trr  be  the  stress  components.  The 
Laplace  transform  with  respect  to  time  of  these  quantities  is  denoted 
by  an  overbar.  Consider  the  problem  where  the  pressure  point-source 
emits  a  step  function  of  strength  M0  at  time  t  ae  0,  where  Mq  =  PQh  and 
where  PQ  Is  the  pressure  amplitude  that  would  occur  in  an  unbounded 
fluid  at  distance  h  from  the  point-source.  Then  the  integral  expressions 
for  the  transforms  of  the  notion  at  r,  z  are 
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where  S(p)  R(p)  +  //J.  »  and  R<P)  =  ('’Z,  “  P**)2  ~  4'^s<tp  P2* 

R(p)  le  the  characteristic  Rayleigh  denominator  in  Lamb’s  problem,  and 
S(p)  is  the  characteristic  Stoneley  (or  Scholte)  denominator  in  Stride's 
problem.  The  path  of  integration,  P  ,  runs  parallel  to  the  imaginary 
axis  between  p  =  0  and  the  smaller  of  p  =  v-*  and  p  =  vf_1. 

By  considering  a  step  function  source,  integrals  have  been  obtained 
in  which  Laplace  transform  variable,  s,  appears  only  in  the  exponentials 
and  the  argument  of  the  Bessel  function.  In  this  case  the  transformation 
to  the  time  domain  is  easily  made  by  operational  means.  The  technique 
is  illustrated  by  a  simple  example. 

Consider  an  integral  of  the  form 

Xfsi  *  F(?)  C  T‘  <2) 


where  F(p)  is  the  ratio  of  polynomials  in  p,  ^ .  7,, 
function  Kn(-spr)  is  the  transform  of 


and 
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is  the  transform  of 
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Thus  the  Integral  1(a)  la  the  transform  of 

j(t)m  _L  f  1*Vrrj) 

Tl,Jn  [(t~  p'r'J  >/lm  <3 


Let  X  =  -pr  +  ^  Bz  ♦  '’l  11  and  d«^orm  the  path,  P  ,  In  the  left  half-plane 
such  that  on  the  del orated  path  T  la  real.  The  deformed  path  la  aynaetrlc 
with  respect  to  the  real  axis  as  In  the  real  part  of  the  Integrand  so 
the  expression  for  l(t)  can  be  written 

l(t>-  TUfc'P(f)  "  '  «> 

where  p+  Is  that  part  of  the  deformed  path  that  lies  In  the  second 
quadrant . 


Changing  the  variable  of  Integration  to  T  gives 


where  p  Is  now  regarded  as  a  function  of  Z  ,  and  ZQ  la  that  value  of  X 
for  which  the  path  of  Integration  touches  the  real  p  axis. 

a 

The  problem  then  is  reduced  to  the  numerical  evaluation  of  an  integral 


(6) 
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If  Z0  corresponds  to  a  real  value  of  p  where  there  is  a  branch  cut,  the 


lower  Halt  of  the  path  of  integration  la  extended  toward  the  origin 
to  the  branch  point  neareat  the  origin.  That  value  of  T  at  the  branch 
point  ia  denoted  by  Tc ,  and 


r.  / 

I/'t)  *  2.  fctx  61-tj  Z)  +  £  [etzGft,  z) 
TT  Jrc  n  Jr. 


(7) 


If  t  <  ln  (6)  or  t  <  7C  in  (7),  I(t)  =  0. 

For  either  (6)  or  (7)  it  la  necesaary  to  evaluate  numerically  one 
or  more  integrals  of  the  type 


■a  af- 

7 r  Jr, 


a/t  6(Z) 


(8) 


The  method  used  ie  the  Gauss- Legendre  method.  Flrat  the  following 
transformation  is  made 


r*  [/  v-  V)^9  *  ] 


from  which  ia  obtained 


T /* 


9  (4*9 


Tl* 


Z/9) 


5  _p J'cteHf) 


(9) 


Next  the  range  of  integration  la  divided  into  K  equal  parts,  /  ;  Zm  N/2K; 


w*-  / 
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To  this  point  no  approximations  have  been  nade. 
In  the  Gauss-Legendre  method  an  integral 


l'  = 

P 


is  approximated  by 


J 


Cll) 


J  */ 

where  Xj  is  the  Jth  root  of  the  Legendre  polynomial  Pj(x), 

and  a 4  la  the  jth  weight. 

J  .  / 


-  1  <  X  <  1, 


_  r  PT 

iJ  ~  Pj<Ki '  J-,  * '  *j 


(12) 


Several  tables  of  weights  and  roots  are  available. 

The  expression  for  the  transformed  Integral  becomes 

I-£Ht  ^ 

j  - » 

The  transformation  procedure  described  here  has  been  applied  to  the  Integral 
expressions  (1.1)  -  (1.5)  and  a  FORTRAN  program  exists  for  the  numerical 
evaluation  of  the  approximating  sums  represented  by  (13). 

The  five  quantities  of  interest  are  presented  in  dimensionless 
fora  as  follows: 

horizontal  velocity,  q  =  /***/*0yt 
vertical  velocity,  w  «  ^uW/P0vf 
vertical  stress,  *zz  “  W»o 
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shear  stress, 

^zr  “  ^zr^o 

horizontal  stress, 

”^rr  -  Tr/>o 

The  following  compilation  defines  terms  used  throughout  the  above 
discussion. 

=  density  of  fluid 
j?*  =  density  of  solid 
V|  =  acoustic  sound  speed  in  fluid 
Vp  =  compressional  speed  in  solid 
vg  =  shear  speed  in  solid 
p  =  a  variable  of  integration 
s  =  Laplace  transform  variable 

=  (vj2  -  p2)l//2  These  expressions  are  multivalued 

functions  of  p.  To  render  them. and 

-ip  -  <»;  -  >2> 

expressions  involving  them,  slnglei 

2 

*1  g  =  (v^  -  p  )  valued,  branch  cuts  are  introduced, 

in  the  case  of  ^  ,  from  p  =  vj*  to 
p  =  «o  and  from  p  =  -v^  to  p  =  -  oo  . 

K  ( o< )  a  modified  Bessel  function  of  the  second  kind,  or  order  n  and 
argument  o(  . 

ytt  =  rigidity  modulus  of  solid. 
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The  results  obtained  above  represent  a  solution  to  the  linear  elastic 
wave  equation.  In  order  to  use  the  results  to  study  the  problem  of  blast 
loading  an  alr/earth  interface,  several  simplifying  assumptions  must  be 
made.  First,  it  was  decided  to  consider  the  response  at  four  ranges, 
corresponding  to  pressure  levels  of  10,000,5000,  1000,  and  500  psi. 

The  parameter  PQ  should  be  larger  than  the  largest  pressure  considered 
and  was  arbitrarily  chosen  to  be  20,000  psi.  That  is  to  say,  at  a 
distance  h  from  the  point  source  the  pressure  would  be  20,000  psi  when 
the  solid  half-space  is  absent.  The  ranges  corresponding  to  the  four 
pressure  levels  considered  were  taken  from  the  data  of  Brode  (RAM),  P-1951, 
March  1960).  Then  for  each  pressure  level  the  parameter  h  was  chosen 
such  that  the  incident  pressure  field  would  have  the  proper  value.  For 
a  pressure  P  at  range  r  we  have 

P  =  hPj/Ch2  4  r2)l/2  t 

For  a  given  P  and  r  and  with  PQ  =  20,000  psi  h  can  be  determined.  Thus 
h  is  different  for  each  pressure  level  and  can  be  regarded  as  the  "effective” 
height  of  burst. 

The  horizontal  blast  wave  velocity  was  taken  from  the  data  of  Brode 
and  from  it  an  "effective"  fluid  speed  can  be  calculated  for  each  range. 

Finally,  the  density  of  the  fluid  is  taken  to  be  the  density  behind 
the  blast  wavefront  in  the  Brode  data  and  therefore  is  also  a  function 
of  range.  This  last  assumption  is  not  very  Important  since  the  fluid/ 
solid  density  contrast  is  always  very  large. 
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On  this  basis  ssvsrsl  problems  were  solved.  Oround  materials  con¬ 
sidered  were  "rock"  (dilatational  velocity  Cj,  =  15,000  ft/sec,  shear 
velocity  ca  =  8,000  ft/sec,  unit  weight  X  a  170  lb/ft3),  "tuff" 

(cp  =  7,000  it/sec,  c8  =  4,000  ft/sec ,  „  130  lb/£t3),  "sediment" 

(cp  =  2,400  ft/sec,  ca  =  800  ft/sec,  =  120  ib/ft3)  and  "dirt" 

(cp  =  800  it/sec,  cg  =  200  ft/sec,  Y  =  100  lb/ft3).  Pressure  levels 
considered  were  10,000,  5,000,  1,000  and  500  psl;  depths  investigated 
were  0,  40  and  80  feet.  Horizontal  and  vertical  normal  stresses,  shear 
stress,  and  horizontal  and  vertical  velocities  were  computed.  Results 
for  tuff  are  included  in  this  presentation;  results  for  rock  and  sediment 
are  given  in  Appendix  A.  Extreme  requirements  of  computer  time  prevented 
extending  the  curves  for  dirt  to  times  at  which  the  Rayleigh  wave 
phenomena  would  occur;  therefore,  no  results  for  dirt  are  given. 

The  Rayleigh  wave  contribution  is  not  apparent  from  examination  of 
the  stress  and  velocity  curves;  the  most  satisfactory  method  of  identifying 
the  Rayleigh  wave  is  by  examination  of  particle  motions.  The  Rayleigh 
wave  is  unique  in  that  particle  motion  is  retrograde  with  respect  to  the 
solid  half-space.  Particle  paths  are  shown  in  Figure  12  for  surface 
points  at  the  four  overpressure  levels.  Since  the  velocity  curves  of 
Figures  10  and  11  begin  at  the  time  corresponding  to  the  peak  stress 
traveling  directly  through  the  fluid  (l.e. ,  at  188  ns)  early- time  motions 
are  not  shown.  The  vertical  velocity  increases  abruptly  after  the  fluid 
pulse  arrival  while  the  horizontal  velocity  decreases  equally  rapidly. 
Positive  directions  of  vertical  and  horizontal  velocities  are  upward  and 
outward,  respectively.  Therefore,  the  motions  following  the  arrival  of 
the  fluid  wave  are  retrograde  and  identify  the  passage  of  the  Rayleigh  wav4. 
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Onca  this  wave  la  located  for  one  pressure  level  It  can  be  Identified 
quite  readily  at  the  other  pressure  levels  as  well;  therefore  It  is 
necessary  only  to  examine  velocity  curves  at  one  overpressure  level. 

For  tuff  the  sharp  downward  thrust  of  the  horizontal  stress  curves 
Immediately  following  the  initial  peak  Is  caused  by  the  Rayleigh  wave. 

The  relative  importance  of  the  Rayleigh  wave  at  the  four  overpressure 
levels  can  be  seen  quite  clearly.  At  500  pel  (nominal)  the  Rayleigh  wave 
horizontal  stress  is  of  much  greater  absolute  magnitude  than  that  from 
the  fluid  pulse.  However,  since  the  stresses  from  these  two  wave  arrivals 
are  opposite  in  sign  the  net  result  of  these  effects  is  a  positive  stress 
arising  from  the  fluid  pulse.  As  overpressures  Increase  the  Rayleigh 
wave  horizontal  stresses  become  less  Important;  for  PBQ  =  5,000  and  10,000 
psl  the  Rayleigh  wave  effect  is  not  of  sufficient  magnitude  to  cause  these 
stresses  to  become  positive. 

The  Rayleigh  wave  has  far  less  effect  on  vertical  stress  than  on 
horizontal  stress.  For  surface  stresses  no  effects  are  shown  whereas 
the  effects  for  depths  of  40  and  80  feet  are  leas  than  10  percent  greater 
for  PRO  *  500  psi  and  less  than  3  percent  greater  for  Pgo  *  10,000  psl. 
Shear  stress  effects  become  quite  large,  particularly  at  lover  pressure 
levels,  as  evidenced  by  Figure  9. 

The  pattern  of  stresses  which  follows  the  Rayleigh  wave  peak  hori¬ 
zontal  stress  is  at  least  partially  obscured  by  oscillations  of  the 
numerical  evaluation.  This  is  particularly  apparent  in  the  horlsontal 
stress  curves  for  Pgo  »  5,000,  1,000  and  500  psl.  Similar  oscillations 
appear  in  the  vertical  and  shear  stress  curves  and  also  in  the  vertical 
and  horizontal  velocity  curves.  Such  fluctuations  could  havs  been 
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eliminated  by  choosing  shorter  time  steps  at  later  times.  This  was  not 
undertaken  In  the  present  study  because  these  oscillations  did  not 
Impair  the  usefulness  of  the  results  In  defining  Rayleigh  wave  effects. 

After  the  passage  of  the  Rayleigh  wave,  vertical  stresses  for 
depths  of  40  and  80  feet  tend  toward  the  values  of  surface  stresses. 
Surface  stresses,  in  turn,  have  stabilized  at  a  value  equal  to  the 
const ant -valued  fluid  pressure  at  these  later  times.  However,  horizontal 
stresses  for  all  depths  Increase  negatively,  by  a  large  amount  and,  In 
the  cases  of  Pso  =>  10,000  and  5,000  psl,  become  larger  than  at  the 
time  of  passage  of  the  fluid  pulse.  For  PBO  =  10,000  psl  these  later 
stresses  are  almost  40  percent  greater  than  those  associated  with  the 
fluid  compress lonal  wave.  It  can  be  established  that  these  late  time 
horizontal  stresses  are  not  a  result  of  the  Rayleigh  wave  disturbance 
by  examining  the  particle  motion  curves  of  Figure  12.  The  transition 
from  retrograde  to  prograde  motion  can  be  seen  clearly.  As  indicated 
in  this  figure  the  change  occurs  at  a  time  corresponding  to  that  of  the 
peak  Rayleigh  wave  stress  for  each  pressure  level  considered. 

For  the  most  part  only  qualitative  conclusions  can  be  drawn  from 
the  above  results  (and  those  given  In  Appendix  A)  because  of  the  simpli¬ 
fications  Involved  in  the  analysis.  The  most  Important  of  the  simplifying 
assumptions  is  that  of  a  non-decaying  pulse.  Were  the  pressure  eource 
a  decaying  pulse  similar  to  the  Brode  air  blast  pulse,  vertical  stresses 
would  exhibit  a  sharp  decay  following  the  initial  peak.  The  horizontal 
stresses  also  would  show  such  a  decay  but  this  would  be  seen  primarily 
at  times  after  the  occurrence  of  the  Rayleigh  wave  peak  stress.  It  is 
unlikely  that  a  decaying  pulse  would  give  rise  to  great  changes  in  the 
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relative  magnitudes  of  the  first  peak  and  the  Rayleigh  wave  peak  stress. 
TOierefore,  the  above  comments  regarding  the  relative  Importance  of  the 
Rayleigh  wave  at  various  overpressure  levels  appear  to  be  valid. 

On  this  basis  It  may  be  stated  that,  for  stresses  In  tuff  at  least, 
the  Rayleigh  wave  Is  of  less  Importance  than  are  the  other  ”alr"-induced 
effects  In  the  close-ln  region  (l.e.  ,  above  1,000  pel).  This  statement 
must  be  qualified  by  adding  that,  taken  alone,  the  Rayleigh  waye  effects 
surpass  the  air- Induced  effects;  however,  the  phasing  of  the  two  effects 
Is  such  that  the  Rayleigh  wave  stresses  are  opposite  In  sign  to  those 
which  exist  prior  to  the  Rayleigh  wave  arrival.  Therefore,  the  net 
effect  of  the  Rayleigh  wave  is  a  reduction  In  stress  rather  than  an  Increase. 
Since  this  phasing  Is  the  dominant  factor  In  determining  Rayleigh  wave 
effects  these  results  are  at  odds  with  the  assumption,  on  which  some 
earlier  analyses  were  based,  that  the  Rayleigh  wave  effects  may  be 
described  independently  of  the  other  alr-lnduced  effects. 
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The  cases  studied  included  four  ground  materials  (rock,  tuff,  sediment 
and  dirt),  four  overpressure  levels,  (10,000,  5,000,  1,000  ^uid  500  pel) 
and  three  depths  (0,  40  and  80  feet).  Properties  assigned  to  the  ground 
materials  were  in  terms  of  dilatatlonal  wave  velocity,  cp,  and  shear 
wave  velocity  cg.  Values  used  were  as  follows:  for  rock,  Cp  a  15,000 
ft/sec,  cg  =  8,000  ft/sec;  for  tuff,  cp  a  7,900  ft/sec,  cg  a  4,000  ft/sec; 
for  sediment,  Cp  =  2,400  ft/sec,  cg  a  800  ft/sec;  for  dirt,  Cp  a  800  ft/sec 
cg  a  200  ft/eec.  Corresponding  values  of  Poisson's  ratio  are  0.30  for 
rock,  G.33  for  tuff,  0.43  for  sediment,  and  0.47  for  dirt.  Results,  in 
the  form  of  time  histories  of  vertical  and  horizontal  stresses  and 
velocities,  are  given  here  for  rock  and  sediment.  Results  for  tuff  are 
given  in  the  main  text.  No  results  are  given  for  dirt;  excessive  com¬ 
putation  time  requirements  prevented  extending  these  results  to  times 
at  which  the  Rayleigh  wave  disturbance  would  occur. 

The  timewise  location  of  the  Rayleigh  wave  is  best  detected  from 
inspection  of  the  particle  motion,  which  is  uniquely  retrograde  during 
passage  of  the  Rayleigh  wavs.  This  has  been  accosqplished  primarily 
from  direct  study  of  velocity  data  although  motion  calculations  were 
helpful.  Identification  of  the  Rayleigh  wave  by  means  of  notion  graphs 
is  discussed  in  the  main  text  and  shown  in  Figure  12. 

Results  for  rock  and  sediment  are  shown  in  Figures  A-l  through  A-10 
and  A-ll  through  A-20,  respectively;  those  for  tuff  are  given  in  the 


main  text 
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Examination  of  the  rock  particle  velocity  curves  of  Figures  A-9 
and  A-10  for  P#0  =  1,000  psl  shove  that  the  retrograde  motion  associated 
with  the  Rayleigh  wave  pulse  begins  at  about  192  ms,  at  which  time  the 
vertical  particle  velocity  changes  sign.  In  Figure  A-7,  in  which  are 
presented  results  for  horixontal  stresses  for  PBO  »  1,000  psl,  the 
pulse  thus  identified  can  be  seen  to  peak  at  200  sui.  At  500  psl 
(Figure  A-8 )  the  peak  can  be  seen  to  occur  at  260  ms,  whereas  the  curves 
for  10,000  and  5,000  psi  (Figures  A-5  and  A-6)  show  the  peak  to  occur 
at  80  and  100  ms,  respectively.  In  these  figures  the  Rayleigh  wave 
contribution  to  the  horiaontal  stress  appears  to  have  a  somewhat  more 
rounded  shape  than  in  Figure  A-7  and  A-8.  This  is  largely  a  matter  of 
the  scales  chosen  for  the  abscissas  of  these  graphs;  drawn  to  the  same 
scale  the  curves  of  Figures  A-6  and  A-7  would  appear  much  less  rounded, 

It  is  apparent  from  these  figures  that  the  relative  Importance  of 
the  Rayleigh  wave  contribution  to  horizontal  stress  in  rock  is  most 
Important  at  the  lower  overpressures.  In  each  case  the  Rayleigh  wave 
produces  horizontal  stresses  opposite  in  direction  to  those  of  the 
pulse  which  is  associated  with  the  passage  of  the  fluid  pulse  front. 

These  horizontal  stress  curves  for  tuff  and  sediment  also  show  decreasing 
Rayleigh  wave  importance  with  increasing  pressure  and  Rayleigh  wave 
stresses  opposite  in  direction  to  those  frosi  the  dilatational  wave. 

It  is  useful  to  discuss  further  the  relstlve  velocities  of  the 
fluid  and  solid  elastic  waves.  For  rock  the  fluid  pulse  is  super seismic 
(i.e. ,  the  point  of  intersection  of  the  fluid  wave  front  ant  the  solid 
surface  travels  at  a  velocity  greater  than  the  propagation  speed  of 
dilatational  waves  in  the  solid)  for  nominal  pressure  levels  of  10,000 
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and  5,000  pal.  For  ?B0  =»  1,000  pal  tha  fluid  pulaa  la  transseismic 
(l.e. ,  lta  velocity  Ilea  between  that  of  dilatatlonal  and  ahear  waves 
In  the  aolid)  and  for  P<0  a  500  pal  It  la  subset sale  (l.e.,  lta  velocity 
la  leaa  than  that  of  ahear  wavea  In  the  aolid).  At  both  1,000  pal  and 
500  pal  the  fluid  pulae  velocity  la  relatively  cloae  to  that  of  the 
ahear  wave  velocity;  the  fluid  and  dilatatlonal  wave  velocitlea  are 
fairly  close  for  PRO  =  10,000,  5,000  and  1,000  pal  and  tranasolaalc 

for  P  _  =■  500  pal.  The  dilatatlonal  and  fluid  wave  velocities  are 

ao 

very  cloaely  matched  at  500  pal  and  relatively  cloae  at  1,000  pal. 

The  fluid  pulae  la  strongly  superseismlc  in  all  cases  for  sediment. 

In  terms  of  the  relative  velocities  it  may  be  stated  from  the 
horizontal  stress  results  that  the  Rayleigh  wave  contribution  la  not  of 
great  importance  In  the  auperselamic  situation  even  when  the  fluid 
velocity  closely  approaches  that  of  dilatatlonal  waves  in  the  solid. 

The  contribution  Is  large  In  both  the  transselsmlc  and  subseismic 
regions  and  reaches  lta  greatest  importance  when  the  fluid  wave  velocity 
closely  approaches  that  of  the  solid  shear  wave. 

Horizontal  stress  curves  for  PBO  *  10,000  and  5,000  pal  for  both 
rock  and  tuff  show  a  rapid  rise  to  a  first  peak*  associated  with  the 
passage  of  the  fluid  pulse.  This  Is  followed  by  a  steep  decay  upon 
arrival  of  the  Rayleigh  wave  followed  In  turn  by  a  later  increase  as  the 

e 

Rayleigh  wave  passes.  The  horizontal  stress  curves  for  tuff  at  1,000 
and  500  psl  also  follow  this  pattern.  The  first  arrival  shown  in 
Figure  A-7,  for  rock  at  1,000  pal  is  not  that  of  the  fluid  equivalent 
of  the  air-blast  shock  front.  The  arrival  time  of  this  early  distur¬ 
bance  corresponds  to  that  of  a  dilatatlonal  wave  Induced  directly 
*In  this  discussion  a  more-negative  (l.e.,  compressive)  stress  is  con¬ 
sidered  greater  than  a  less-negative  streas. 
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beneath  the  burat.  The  fluid  wave  front  arrives  at  188  ns,  when  the 
pulse  has  gone  positive  by  an  extremely  large  amount.  The  Rayleigh 
wave  peak  occurs  at  200  ms  and  gives  the  most  positive  point  on  the 
curve  for  each  depth.  In  Figure  A-8,  for  rock  at  500  psl  the  fluid  pulse 
arrives  after  the  positive  peak;  the  effect  of  thin  pulse  can  be  seen 
in  the  abrupt  change  of  the  surface  stress  from  4 .350  to  2 .300  psl 
in  about  10  ms.  In  this  case  the  Rayleigh  wave  peak  occurs  at  280  as 
and  the  first  arrival  is  that  of  a  dllatational  wave  Induced  directly 
beneath  the  pulse.  It  la  apparent  that  the  surface  disturbance  has 
outrun  the  fluid  wave. 

The  critical  situation  in  rock  is  shown  in  Figure  A-7  for  P-0  a  1,000 
psl  where  the  fluid  pulse  is  passing  from  the  transseisalc  to  the  subseismic 
case.  This  only  occurs  in  the  region  of  Interest  for  materials  of  high 
seismic  velocity  such  as  the  rock  considered  here;  even  for  tuff  the 
location  of  this  transition  is  far  beyond  the  distance  for  which  Pso  ■ 

1,000  psl.  8uch  high-velocity  materials  also  possess  great  strength 
and  high  modulus  of  elasticity.  Such  materials  are  far  more  able  to 
withstand  the  large  stresses  associated  with  the  Rayleigh  wave  in  the 
transseisalc  region  than  are  softer  materials.  Thus  horizontal  stress 
effects  appear  to  be  largest  in  materials  which  are  beat  suited  to 
supporting  such  stresses.  There  is,  of  course,  some  question  as  to 
the  significance  relative  to  design  of  the  hlgh-intenslty  tensile 
stresses  such  as  those  of  Figure  A-7  (l.e. ,  4,000  psl). 

The  absolute  magnitude  of  horizontal  stress  in  rock  arising  from 
the  passage  of  the  Rayleigh  wave  is  in  no  Instance  greater  than  the  peak 
fluid  Induced  horizontal  stress  at  PBO  =  10,000  psl.  For  tuff  the  Rayleigh 
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wave  stress  Is  nowhere  larger  than  the  fluid-induced  stress  at  P  * 

so 

1,000  psi.  Results  for  sediment  are  incapable  of  supporting  such  a 
statement  directly  because  the  computations  were  halted  before  the 
Rayleigh  wave  peak  occurred  for  P|o  a  1,000  and  500  psi.  Comparison 
of  the  sediment  curves  for  horizontal  velocity  with  those  for  rock 
and  tuff  at  PSQ  =  10,000  and  5,000  psi  indicates  that  the  Rayleigh 
wave  effects  for  sedlaent  are  less  Important  relative  to  the  fluid- 
induced  effects  than  are  the  effects  in  rock  or  tuff. 

Surface  vertical  stress  curves  for  rock  at  all  pressure  levels 
show  an  abrupt  rise  on  arrival  of  the  fluid  pulse  and  a  constant 
value  of  stress  thereafter.  Stress  curves  for  depths  of  40  and  80  feet 
closely  approximate  the  surface  curves  for  P|0  »  10,000  and  5,000  psi; 
for  Pso  =  1,000  and  300  psi  the  40  and  80  feet  curves  differ  considerably 
from  the  surface  curves  at  earlier  times.  These  differences  at  early 
times  are  attributable  to  the  outrunning  ground  shock  at  these  lower 
pressure  levels.  At  1,000  psi  the  early  stresses  are  Insignificant 
relative  to  the  stresses  which  occur  after  the  fluid  pulse  arrival. 

The  relative  importance  of  the  early  stresses  is  greater  for  P#0  a  500 
psi  but  fluid  pulse  stresses  remain  greater  than  those  associated  with 
the  Rayleigh  wave.  The  Rayleigh  wave  effects  in  each  case  tend  to 
counter  the  stresses  from  the  fluid  pulse.  At  later  tints  curves  for 
all  depths  tend  toward  the  value  of  surface  pressure,  as  is  to  be  expected. 

Curves  for  vertical  stresses  in  tuff  and  sediment  exhibit  the  same 
patterns  as  those  for  rock  with  the  exception  that  no  outrunning  occurs 
for  these  softer  materials.  Sediment  results  for  later  times  at  9mo  » 
1,000  and  500  psi  are  unavailable  but  the  results  appear  to  follow  those 


for  tuff 
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Velocity  curves  for  rock  are  shown  In  Figures  A-9  and  A-10  for 
vertical  and  horizontal  velocity,  respectively.  For  sediment  the 
corresponding  figures  are  A-19  and  A-20  and  for  tuff,  Figures  10  and 
11  of  the  main  text.  Motion  phenomena  are  strongly  affected  by  the 
shape  of  the  pulse  and  no  coonaents  regarding  magnitudes  of  the  velocities 
are  in  order.  These  curves  are  useful  primarily  in  identifying  the 
Rayleigh  wave  from  the  retrograde  motion.  This  motion  can  be  found 
readily  by  considering  simultaneously  the  changes  in  vertical  and 
horizontal  motion. 

Because  of  the  simplifications  involved  in  this  analysis  of  the 
Rayleigh  wave  only  limited  conclusions  are  drawn.  These  are  sufficient 
to  define  the  relative  importance  of  the  Rayleigh  wave  contribution, 
at  least  with  regard  to  stresses.  These  conclusions  are  as  follows: 

1.)  The  Rayleigh  wave  is  never  of  primary  importance  for  horizontal 
stresses  in  the  superseismlc  case;  therefore,  only  materials  having 
seismic  velocities  greater  than  8,300  ft/sec  (the  air-blast  wave  velocity 
at  PBO  =  1,000  psl)  will  exhibit  important  Rayleigh  wave  effects  in 
the  close-in  region.  2.)  The  Rayleigh  wave  is  never  of  great  importance 
for  vertical  stresses  at  any  location  for  any  material,  3.)  The 
Rayleigh  wave  horizontal  stresses  are  opposite  In  sign  to  those  accompanying 
the  air  blast  wave  and  tend  to  reduce  the  effects  of  the  latter;  for 
vertical  stresses  the  results  generally  are  additive  but  the  Rayleigh 
wave  contribution  Is  of  minor  importance.  4.)  When  Rayleigh  wave 
effects  are  important  the  stresses  involved  are  tensile. 

Based  on  the  above  conclusions  stresses  arising  from  the  Rayleigh 
wave  appear  to  be  of  less  importance  than  are  other  air- induced  effects 
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in  the  close-in  region  and  can  be  neglected  tor  Boat  design  application* 
without  introducing  appreciable  error  into  the  analysis.  In  order  to 
make  a  comparable  statement  regarding  ground  motions  it  will  be  necessary 
to  base  any  future  analysis  on  a  more  realistic  representation  of  the 
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Vertical  Stress  In  Tuff 

Pso  «*  1,000  psl,  Seismic  Velocity  =  7,900  ft/sec 
Shear  Wave  Velocity  s  4,000  ft/sec. 
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Vertical  Stress  in  Tuff 

Pso  =  500  psi,  Seismic  Volocity  =  7,900  ft/soc 
Shear  Wave  Volocity  =  4,000  ft/sec 
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Horizontal  Stress  in  Tuff 

Pao  *  9,000  psl,  Seismic  Velocity 

Shear  Wave  Velocity  *  4,000  ft/sec 
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Vertical  Stress  In  Rock 

Pso  3  300  psi*  Soismlc  Velocity  =  15,000  ft/soc 
Shear  Wave  Velocity  =  8,000  ft/soc 
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Vertical  Particle  Volocity  for  Rock 

Pgo  3  1,000  psl,  Seismic  Velocity  =  15,000  ft/sec 

Shear  Wave  Volocity  3  8,000  ft/sec 
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FIGURE  A- 11 

Vortical  Stress  In  Sediment 

P3Q  =»  10,000  psi  ,  Seismic  Velocity  =  2,400  ft/sei 
Shear  Wave  Velocity  =  800  ft/soc 
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FIGURE  A-18 
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Horizontal  Stress  in  Sediment 

pso  =  500  pal,  Seismic  Velocity  =  2,400  ft/soc 
Shear  Wave  Velocity  =.  800  ft/sec 
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Horizontal  Particle  Velocity  for  Sedlnent 

Pjo  ■  1,000  pal,  Selaalc  Velocity  ■  2,400  ft/eec 

Shear  Wave  Velocity  a  800  ft/eec 
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